Intrinsic Spin Hall Effect in the Two Dimensional Hole Gas 
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We show that two types of spin-orbit coupling in the 2 dimensional hole gas (2DHG), with and without inversion 
symmetry breaking, contribute to the intrinsic spin Hall effect [J, |2|. Furthermore, the vertex correction due to impurity 
scattering vanishes in both cases, in sharp contrast to the case of usual Rashba coupling in the electron band. Recently, 
the spin Hall effect in a hole doped GaAs semiconductor has been observed experimentally by Wunderlich et al 0. From 
the fact that the life time broadening is smaller than the spin splitting, and the fact impurity vertex corrections vanish 
in this system, we argue that the observed spin Hall effect should be in the intrinsic regime. 

PACS numbers: 73.43.-f,72.25.Dc,72.25.Hg,85.75.-d 



Recent theoretical work predicts dissipationless spin 
currents induced by an electric field in semiconductors 
with spin-orbit coupling 0, 0, Q. The spin current is re- 
lated to the electric field by the response equation 



J i 



E k 



(1) 



where fjj is the current of the i-th component of the spin 
along the direction j and e^-fc is the totally antisymmet- 
ric tensor in three dimensions. Because both the electric 
field and the spin current are even under time reversal, 
the spin current could be dissipationless or intrinsic, in- 
dependent of the scattering rates. The response equa- 
tion (JTJ was derived by Murakami, Nagaosa and Zhang Q 
for p-doped semiconductors described by the Luttinger 
model of the spin-3/2 valence band. In another proposal 
by Sinova et. al. |2|, the spin current is induced by a 
in-plane electric field in the 2-dimcnsional electron gas 
(2DEG) described by the Rashba model0. 

The spin Hall effect predicted by these recent theoret- 
ical works is fundamentally different from the extrinsic 
spin Hall 0, 0] effect due to Mott type of skew scatter- 
ing by impurities. The intrinsic spin Hall effect arises 
from the spin-orbit coupling of the host semiconductor 
band, and has a finite value in the absence of impuri- 
ties. On the other hand, the extrinsic spin Hall effect 
arises purely from the spin-orbit coupling to the impu- 
rity atoms, and it is not a bulk effect like the ordinary 
Hall effect. Because the extrinsic arises only from the 
impurities rather than the host atoms, its magnitude is 
typically many orders of magnitude smaller. The issue 
of impurity contributions to the spin Hall effect has been 
intensively investigated theoretically. Remarkably, Inoue 
et. al. 7] calculated the vertex corrections due to impu- 
rity scattering in the Rashba model of the electron band, 
and found that the vertex correction completely cancels 
the spin Hall effect. Other analytical works have obtained 
similar results || On the other hand, a number of numer- 
ical calculations have shown that the spin Hall effect is 
independent of the disorder in the weak disorder limit 
I& [lO( ■ Currently, this disagreement between the analyt- 
ical and the numerical results is still not settled. In a 
insightful paper, Murakami fllj showed that the problem 
of the vertex correction does not occur in the Luttinger 



model of the hole band 11]. In fact, the vertex correc- 
tion is identically zero, rendering the original prediction 
of Ref. exact in the clean limit. 

Experimental observation of the spin Hall effect has 
been recently reported by Kato et. al m a electron 
doped sample and by Wunderlich et al in a 2DHG[3j. In 
this paper, we analyze the 2DHG experiment. In order to 
firmly establish the intrinsic spin Hall effect, one needs 
to establish two things. First of all, the experimental 
system needs to be in the clean limit, which is the case 
of the 2DHG experiment, as shown in the experimental 
paper [3- Secondly, one needs to show that the effect is 
robust in the clean limit, not cancelled by the vertex cor- 
rections due to impurities. We shall show that the spin 
Hall effect in the 2DHG arises from two contributions, 
one from the Luttinger Hamiltonian describing the split- 
ting between the light and the heavy hole bands, and one 
from the structural inversion symmetry breaking (SIA) 
of the 2DHG band, with a spin splitting scaling as fc 3 
[I3I IhH . The later form the the spin-orbit coupling has 
been studied by Schliemann and Loss 14] in connection 
to the intrinsic spin Hall effect. This is different from 
the Rashba Hamiltonian of the 2DEG band, where the 
spin splitting scales with k. Remarkably, we find that 
the vertex correction due to impurity scattering vanishes 
for both types of spin-orbit couplings in the 2DHG band, 
in sharp contrast to the case of 2DEG. While the cal- 
culation details are complicated, the intuitive reason is 
simple: The two types of current vertices in the 2DHG 
have p wave and d symmetries, respectively. When these 
current vertices are averaged over the s wave impurity 
scatters, the vertex corrections vanish. These two key 
facts establish a firm foundation to interpret the recent 
experiment by Wunderlich et al in terms of the intrin- 
sic spin Hall effect, where impurities play an unessential 
role. 

The Hamiltonian for a 2-dimensional hole gas is a sum 
of both Luttinger and spin-5 = 3/2 SIA terms: 



H=( ll + ^ 2 )^~^(k-S) 2 + a(Sxk)-z (2) 
2 2m m 

where the confinement of the well in the z direction makes 



2 



the momentum be quantized on this axis. The crucial 
difference between the SIA term for 2DHG and Rashba 
term for the 2DEG lies in the fact that S is a spin 3/2 
matrix, describing both the light (LH) and the heavy 



(HH) holes. For the first heavy and light hole bands, 
the confinement in a well of thickness a is approximated 
by the relation < k z >= 0, (k 2 ) s» (nh/a) 2 . The energy 
eigenstates are: 



E 



hh _ 7i , 2 , 1 . 



2m 



k 2 ± ^-ak - \la 2 k 2 ± ^-k(k 2 + (k 2 )) + -%(fc 4 + (k 2 ) 2 - k 2 (k 2 )) 

E LH = 7l_ fc 2 ± l ak + I 2 k 2 ± ?Hl k ( k 2 + / fe 2\) + jL( fc 4 + (fc|\2 _ fc2(jfc2\) 

2m 2 V m m 2 



(3) 



The heavy and light hole bands are split at the T point by 
A = 272 (fc 2 ) /m [la llfij. Depending on the confinement 
scale a the Luttinger term is dominant for a not too small, 
while the SIA term becomes dominant for infinitely thin 
wells, which correspond to high junction fields. 

By expanding the above formulas for small k << (k z ) 
it is seen that the spin splitting of the HH bands is fc 3 
whereas the spin splitting of the LH bands is fc, in agree- 
ment with dl G3 
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(4) 
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Figure [1] gives a typical band structure for GaAs (71 = 
6.92, 72 = 2.1) with a T point gap of AQmeV and a Fermi 
momentum splitting of the hole band at Fermi momen- 
tum (0.2n.m ) of hmeV , which require a SIA splitting 
a 10 5 m/s . 

We can expand the second term in the anisotropic Lut- 
tinger hamiltonian in terms of Clifford algebra of Dirac 
T matrices {r a ,r 6 } = 2S ab hx4 (a, b = 1, 5)0. Since 
\J (k 2 ) — and (fc 2 ) ^Owe see that the effect of confine- 
ment renders the d a, s of |4|: 



(fc • S) 2 = d a T a - di = 0, d 2 = 0, d 3 
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dA = - — {K - d 5 = --(2(kj) -hi- fc 2 ) 

(5) 

Since calculation with the full Hamiltonian J5J is an- 
alytically impossible, we concentrate on different situ- 
ations which maintain analytic predictability. We now 
consider the case of of small junction field and neglect 
the SIA term: 



H = ^(k 2 x + k 2 + (k 2 z )) + ^d a r a 



72 



2m 



rn 

'„2 1 /u2\ 



(0) 



The energies arc E LHMH = ^{k 2 + (fc 2 )) ± d, (d = 

= v /fc4 + ( k z) 2 - (k 2 z )k 2 ). At fc = the heavy 
and light hole bands are split by a gap of AE — 



FIG. 1: Approximate band structure of the 2DHG (A = 
40meV, and the spin splitting of the heavy hole band at the 
k,F is around 5meV). The confinement produces a F point 
gap between the light and heavy hole bands, whereas the SIA 
produces splitting in the previously degenerate Kramers dou- 
blets - the heavy and light hole bands. 



2^(fc 2 ) w 2^(7rh/a) 2 . In the experiment recently re- 
ported 0, this energy gap is of order AE — AQmeV, 
which corresponds to an a = 8.3nm thick quantum well. 
The value quoted in the experiment is roughly 3 — 4nm, 
making our simplistic prediction rather accurate. 

We would like to compute the response of spin current 
J\ = ^{S l , ||r-} to an electric current Jj = 



Q\Aw n 
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spin conductance 
+0 zjz and S lves: 



(TJl{u)J,)e—du 



is then defined as en 
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where rf b is a tensor denned in n± are the Fermi 
functions of the two bands and e = 2m[k^. + ky + (fcf)) 
is the kinetic energy. The last term is the conserved spin 
conductance (which represents the response of the spin 
projected onto the HH and LH bands |^|), whereas the 



first term is the contribution of the non-conserved part 
of the spin. Upon momentum integration, since motion 
along the z axis is prohibited due to confinement, the 
only non-zero components are o\ 2 = — °2i which yield: 



'12 
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3 7i 
2 2 72 



2(fc 2 + (fc 2 )) 



3^ + (fc 2 ) 2 - (fc 2 )fc 2 



- ln[2\A 4 + (fc 2 ) 2 - {k 2 )k 2 + 2k 2 - (fc 2 )] 



2(fc 2 2 



V^ 4 + (fc 2 ) 2 - (fc 2 )fc 2 



k—k L 



(?) 



i 3(noncons) 

t 3(cons) 

spin conductance c 12 



where k^n , fcffff arc the fcrmi momenta of the light and 
heavy hole bands. For the experimental data 0, the light 
hole band is fully occupied, so k^H = while y/ (fc 2 ) = 
3.7 x 10" 26 fcg m/s and k H H = 3 x 10" 26 fcg m/s. The 
first two terms are due to the non-conserved spin and for 
2^. The last term is the conserved 
= 0.6 x -L. The total spin 
conductance is therefore crf 2 = I^Fj m g°°d agreement 
with the numerical estimate in Q ■ In the case of infinite 
confinement, yj (fc 2 ) — > oo the spin conductance from the 
Luttinger term vanishes, as it should, since we would 
then enter a SI A dominated regime. 

We now investigate the effect of disorder on the Lut- 
tinger spin Hall conductance. In particular, we want to 
find out the vertex correction. The free Green's function 
in our system is defined as Go(k, E) — [E — if] -1 : 



that the free vertex cancels 



G (k,iu n ) 



iui n - e(k) + ^dgYg 
{iuj n - e(k)) 2 - "i\d 2 /n. 



(10) 



We model the disorder as randomly distributed, spin- 
independent identical defects V(r) = uJ2i $( r ~ in 
the Born approximation, the self-energy is related to the 
free Green function Y,{iuj n ) — ni mp u 2 J j4jrp GqQs., iuj n ) 
Since J dkd-j — J <ikc?4 = 0, the self energy is an isotropic 
function of fc: 



E(iw n ) — Tii m pi 



dk 



ioj n - e(k) 
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(2^) 2 (^„-e(k)) 2 - 7 2 d 2 /™ 2 

(11) 

where d§ — — i(2(fc 2 ) — fc 2 ). This is different from the 
bulk Luttinger case, where the ds(fc) integral over fc van- 
ishes as well, but the difference is not essential. The 
full impurity Green function is G(k, iu> n ) = Go(k, iui n + 
E(iu> n )). The current vertex satisfies a Bethe-Salpeter 
equation similar to ^7|. Similar to the case of [Tlj . 
since the Green's function is an even function of the in- 
plane total momentum, while the charge current opera- 
tor Vj = dH/dkj is momentum-odd in the components 
kj (because the Hamiltonian H is even in fc, it turns out 



dk 



G(k, iuj n )Vj{k)G(k, iw r , 



(12) 



And hence the vertex correction which is an iterative 
function of the free vertex vanishes as well ^lj • 

To see the effect of a very small SIA splitting (akp << 
A, akp << h/r) on the Luttinger spin Hall conduc- 
tance, we treat it in perturbation theory. The calculation 
is long and not particularly revealing, so we just give the 
result. Due to the fact that the Luttinger current op- 
erator is odd in fc while the SIA current operator is a 
constant matrix, the first order contribution in a van- 
ishes. The first nonzero contribution is of order a 2 . 

We now turn to the opposite case, of strongly confined 
quantum wells, in which the SIA term is likely to dom- 
inate. We model the system by a Gamma point gap A 
plus a spin 3/2 Rashba term a(k x S)z 0. We compute 
the spin conductance and expand it in terms of the ra- 
tio between the SIA spin splitting and the T point gap, 
< 1. The spin conductance gets a contribution from 
the HH band: 



3(HH) 
T 12 



_9_ n a 2 m HH 
8tt [ 2A ' 



(13) 



For infinitely thin quantum wells, A — > oo, the HH spin 
conductance is 9/87T which is (besides a re-scaling fac- 
tor of 2 in th e sp in current definition), the same as that 
obtained in who studied the hole gas by starting 
with the effective HH Hamiltonian directly. The second 
term in Eq^2] is t ne first order finite thickness correc- 
tion. If the Fermi level is low enough, there is also a 
light-hole band contribution to the spin conductance of 



order a 



3(LH) 



'12 - ^f(1 + 3a 2A H )- The vertex correction 
for a spin-3/2 Rashba- like system has been computed in 
3D perturbatively to first order in a and was found to be 
finite ■ We now compute it exactly for the 2D case in 
the heavy hole band. 

Since working with spin 3/2 matrices is cumbersome 
and we do not need the LH states as they are fully filled 
0, we now project the system onto the heavy hole states 



4 



and work with the truncated Hamiltonian 0, ^| : 



H 



2m 



(14) 



which becomes exact in the limit of very confined quan- 
tum wells. The spin Hall conductance in the disorder- 
free case is 9/87T, as previously obtained in Q and as 
obtained above as the limit of strongly confined quantum 
wells. The Hamiltonian can also be expressed: 



H = h \i(k)ai, i = x,y 

2m 



(15) 



Where Ai = l3k y (3k 2 - k 2 y ) and A 2 
(3 a constant. Let A(fc) = \/\\. 



= /3fcx(3A;2 

The Fermi sphere is 

■ 2 

isotropic since the energy levels are E± = j— ± A. The 
disorder-free Green function is: 



G k, iu) n ) = - > — 

2 ^-f iuj n - E s 

s — ± 



where Ai = A^/A. The self energy for s-wave scatter- 
ing of electrons becomes a state-independent constant 
(not a matrix) £(zu;„) = rii m ' ' k 



/ (^F G o(k, iu n ) where 



riimp is the density of impurities while u is the impurity 
potential strength . Since Ai are odd functions of the 
momentum components fej, the integral J dk/(fc)Ai = 
where f(k) is any isotropic function of k. Since the 
spin orbit coupling small (much smaller than the Fermi 
energy), the density of states at zero order is a con- 
stant D = m/2irh 2 while the ak 3 term in the Hamil- 
tonian contributes to with only a first order correction. 
The full Green function in the presence of impurities is 
G(k, iuJn) = Go{k,iuj n + T,(iuj n )). The spin dependent 
part of the charge current operator Vj(k) — dH/dkj 
turns out to have d-wave symmetry (for example, the 
spin dependent part of the V x operator reads 6(3k x k y a x + 
k x ), with 3/3(fc 2 — k 2 x )a y ) and it vanishes when integrated over the 
isotropic Fermi surface. This is the deep intuitive reason 
as to why the vertex correction cancels in this case, as 
we rigorously show below. By contrast, in the electron- 
band Rashba case, the spin-dependent part of the charge 
operator is a constant. The current vertex function 
Kj(k,iw n , iv m ) = (GQs.,iu}„)Vj(k)G(k,iLJ n + iv m )) is a 
matrix function that does not commute with either the 
charge current operator or the Green's function. (...) is 
an impurity average. It satisfies the Bethe-Salpeter equa- 
tion: 



(16) 



Kj(k, iu n ,iv m ) = G(k,iu n )Vj(k)G(k,iu n - iu m ) + n imp u 2 G(k, iuj n ) j ^ Kj (q, iw n , iv m )G(k, iu n - iv m ) (17) 



Integrating both the right and the left hand side over 
the momentum k, we see that the vertex correction 
AVj(iuj n ,ivm) = J j0^K 3 ((i,iuj n ,w m ) satisfies: 



2 f dk 

+n mp u J 



G(k, iu n )Vj(k)G(k, iw n - iv m ) + 
G(k, iu) n )AVjG(k, iuj n - iv m ) (18) 



Since the vertex correction AVj(iu> n , iv m ) is a 2 x 2 ma- 
trix, it can be decomposed in the basis of the identity 
matrix and the 3 pauli matrices: 



AVj{iu n) iu m ) = 2J A^(iw„, iv m )a^ , n = 0,1,2,3 

(19) 

where er° = ^2x2, the identity matrix, and cr 1 ' 2,3 are the 
3 Pauli matrices. The Aj(iu) n , iv m ) are scalars. By intro- 
ducing the decomposition in the vertex equation, multi- 
plying to the left of both sides of the equal by a <r a matrix 



and taking the trace of the above equation, we obtain: 



2A£ = A?(iw n , iu m ) + A^M^(to„, iv m ) 
r die 

M v » = n imp u 2 J _—Tr[^G(k, uj n )o»G(k, iu n - iv m )\ 

/dk 
—-^Tr[a u G(k, iuo n )V, (k)G(k, iu n - iv m )\ (20) 



By expanding and evaluating (this uses the obser- 

vation that J dkX l (k)X j (k) — 6^ as well as GfGf = 
=j£S(e — E s ), where R,A stand for the retarded and 
advanced Green's functions, and r = fi 3 /ni mp u 2 m) we 
observe that it is diagonal in v, that is M v ^ = b vii . 
Expanding the traces in Eg. 1)20(1 . and since Ao(fc) = 
Aa(fc) = it is easy to observe that (after azimuthal 
integration) A®(iuj n ,iv m ) = A?(iu n , iv m ) = and hence 
the vertex corrections Aj(iui n ,ip m ) = A.j(iuj n ,iv m ) = 0. 
We now have for the vertex correction A^(iu) n ^iv m ) = 
Aj(iu> n , iv m ), v = 1, 2 and j — x,y where after expand- 



5 



ing the traces: 



E 



dk 



27T) 



2(z - E a ){zT - E, 



(21) 

with z = iu; n + E(io; n ), z' = iw„ — iv m + Y,(iuj n — iv m )- 
We now compute this for v = 1, the case = 2 being 
identical. Let j — 1 and we find for the numerator of the 
integrand in EqJ2IJ: 



fc3 



6ss'(3k +6f3k x k y (l -ss') 
(22) 

Upon integration over dk the above expression vanishes 
due to the integral over the azimuthal angle and hence 
A\ {iui n , ifm) = 0. For the case v = l,j = 2, the numer- 
ator of EqJ2U gives: 



life time broadening is smaller than the spin splitting, 
we hence conclude that the spin Hall effect observed in 
-gjTj- 3] should be in the intrinsic regime. The dissipationless 
"spin Hall conductance can be systematically determined 
by extrapolating the ratio of life time broadening and the 
spin splitting to zero. 
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s')—+6ss'pk 
m 



+ 30(k 2 x -k 2 y )(l-ss') 
(23) 

which also vanishes upon azimuthal angle integration 
A\{iuj n ,iv m ) = 0. In an identical way all the compo- 
nents of the vertex correction tensor vanish. 

We have analyzed the spin Hall transport in the case 
of a two-dimensional hole gas. We showed that for rel- 
ative weak confinement the spin-Hall conductance is of 
Luttinger type and is equal to roughly 1.9e/8n for the 
value of parameters in [3j. For strongly confined quan- 
tum wells, the system is dominated by a structural in- 
version asymmetry term of spin-3/2 SlA-type. The spin 
conductance for this system is 9e/87r plus a correction 
dependent on the quantum-well size. We perform the 
full vertex correction and show that it vanishes for both 
Luttinger and SIA cases. This is in striking contrast to 
the fc-linear Rashba case, where the vertex correction is 
of the same magnitude and of opposite sign to the spin 
orbit coupling strength. Coupled with the fact that the 
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